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Abstract 

' Recently developed analytic methods in the framework of the Field Cor- 

relator Method are reviewed in this series of four lectures and results of cal- 
culations are compared to lattice data and experiment. Recent lattice data 
demonstrating the Casimir scaling of static quark interaction strongly sup- 
port the FCM and leave very little space for all other theoretical models, 
e.g. instatnton gas/liquid model. Results of calculations for mesons, baryons, 
quark-gluon plasma and phase transition temperature demonstrate that new 
analytic methods are a powerful tool of nonperturbative QCD along with 
q^. lattice simulations. 
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1 General QCD formalism for Green's functions. 

1. Starting to work in QCD one should have in mind some specific features, which 
make this theory unique; two of these are most important. 

In QCD (as in any local gauge theory in the confining phase) it is impossible 
to introduce asymptotic states of quarks and gluons, hence the popular formalism 
of LSZ ]1|] is not applicable. Therefore one works with off-shell quantities, like the 
Green's functions and later on one introduces Hamiltonian and scattering ampli- 
tudes. 

2. In QCD the Nonperturbative (NP) fields are Euclidean (c.f. ||) and therefore 
one should work out the formalism of Euclidean vacuum averaging first, and later 
one goes over to the Minkowskian space (inevitable for scattering amplitudes.). 

In this section we introduce first simple general and exact formulas for Green's 
functions, Hamiltonians etc. and later discuss possible approximations. 

Consider a general partition function of QCD 

Z{J^ v ,r])= J DA^D^ + e- s(A '^ +)+ I^ A ^ + ^ + ^ d4x (1) 
where S is the standard QCD Euclidean action. 

S=^J trF* u d 4 x - J i) + {id + im + gA)d 4 x. (2) 
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Here gauge fixing terms and Faddeev-Popov ghosts are suppressed for simplicity. 

A spicific feature of local gauge theories is that only gauge-invariants can enter 
as physical states, and local gauge invariance cannot be spontaneously violated (the 
Elitzur theorem). Therefore the only way to proceed is to form gauge-invariant 
Green's functions from ([l]). E.g. for the qq state one might have a local or nonlocal 
meson Green's function. 

G%l(xx\yy) = (x)T ^(x)^ (y)T , or^+I/W+IW) (3) 

where r/,r in = 1, 75, 7^, $(x, y) = P expig J* A^dz^ and (...) means averaging 
over DA IM DipDip + as in (0). 

Integrating over DipDif) + in ([[]) one obtains for the meson Green's function. 

G (M) = / flV" S ' (A) det(m + D){tr a [T f S(x,y)T in S(y,x)}- 

- tr a (S(x,x)Tf)tr a (S(y,y)T in )}, (4) 
with notatins for the quark Green's function 

S(x,y) = (x\(m + D)- 1 \y). (5) 

In QED one can define for G^ the Bethe-Salpeter (BS) or similar equations to take 
into account the proper subclass of perturbative diagrams. In QCD this way is im- 
possible since both BS and Dyson-Schwinger eqs are not gauge-invariant, whenever 
one approximates the kernels by finite sets of diagrams. 

For gluons one should introduce gauge-invariant initial and final states made of 
Ffiv, Dp. E.g. for local states */(» n ) = trF flu (x)F a / 3 (x) or trD p F IIU (x)F a p(x)) etc. 
and define Green's functions 

GfS = (*;(x)* to (y)>. (6) 
Integration over DipDip + allows to obtain only the form 

Gfjl = J DA^e- 8 ^ det(m + D)V+(x)* in {y). (7) 

To proceed one should use the Background Field Formalism (BFF) and the 
Feynman-Schwinger (or world-line) Representation (FSR) which will be discussed 
in the next lectures. 

Here we shall give connection to the Hamiltonian technique, which is useful both 
for bound states and for scattering. 

The following general relation holds true for any Green's function 

G f}n( x i y\ x 'i V') = (/> x > y\e~ HT \in, x', y') = (/, x, y\n > e~ MnT < n\in, x', y') (8) 

with M n , \n > being eigenvalues and eigenf unctions of H respectively. Here we have 
defined a set of hypersurfaces {J2t} 011 which the Hamiltonian is constructed, and 
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the evolution parameter t, perpendicular to J2n which varies between, say, t = for 
the in-state and t = T for the final state. 

The Hamiltonian H is the general QCD Hamiltonian which is known in several 
gauges, e.g. Coulomb 0, temporal (4| or in polar representation [Q. It is a matrix, 
containing creation and destruction of quarks and gluons, and angular brackets in 
(^D do single out the initial and final indices of this matrix. 

The most important feature of QCD is confinement, which means that some part 
of gluons are condensed in the form of strings between charges. Unfortunately the 
general QCD Hamiltonians [f§-|| are not yet helpful in understanding confinement 
properties of QCD, and one cannot use it as a tool in exploring vacuum properties. 

Therefore our main aim would be to formulate effective Hamiltonian for the 
valence (noncondensed) gluons and quarks. 

The way to do it will be described in next lectures. 

2 Background Field Formalism 

One of basic mysteries of QCD is the fact, that gluon field plays two different roles: 

i) gluons are propagating, and at small distances this process can be described 
perturbatively, leading in particular to color Coulomb interaction between quarks 
(antiquarks); 

ii) gluons form a kind of condensate, which serves as a background for the prop- 
agating perturbative gluons and quarks. This background is Euclidean and ensures 
phenomena of confinement and CSB. 

Correspondingly we shall separate the total gluonic field into perturbative 
part a M and nonperturbative (NP) background B^. 

A^ = B^ + a M (9) 

There are many questions about this separation, which may be answered now 
only partially, e.g. what exactly is the criterion of separation. Possible answer 
is that perturbative fields a M get their dimension from distance (momentum), and 
therefore all correlators of fields a M (in absence of By) are singular and made of in- 
verse powers of (x — y) and logarithms, where enters the only dimensional parameter 
of perturbative QCD - Aqcd- Therefore evidently any dimensionful constant, like 
hadronic masses or string tension cannot be obtained as a perturbation series. In 
contrast to that, NP fields B^ have mass dimension due to the violation of scale 
invariance which is intrinsically present in the gluodynamics Lagrangian. The ori- 
gin of separation (|9])) is clearly seen in the solutions of nonlinear equations for field 
correlators || : a perturbative solution of those leads to singular power-like field cor- 
relator, whereas at large distances there is a selfconsistent solution of the equations, 
decaying exponentially with distance with arbitrary mass scale, since equations in || 
are scale-invariant. Full solution including intermediate distances produces mixed 



3 



perturbative-nonperturbative terms, containing both inverse powers of distance and 
exponentials. For these terms criterion of separation fails. 

One can avoid formally the question of separation principle (and of double count- 
ing) using t'Hooft identity |7|], which allows to integrate in (1) independently over 
Bfj, and a a : 

Z = ^-,j DBtf(B)D1>D$Da ll e L ( B+a '> (10) 

Here L(A) = L + L\ + L int , and A u is taken to be B u + a u ; the weight rj(B) is 
arbitrary and may be taken as constant. For the exact formalism starting from 
fllOD we refer the reader to @, §, and here we only quote the form of background 
propagator G^ v of the gluon, which is found from L(B +a) and will be used in what 
follows, 

G au = {D\5 au + 2igF au )- 1 (11) 

where D\ = d\ — igB x , and = d u B u - d u B a - ig[B^ B v \. 

To define the perturbation theory series in ga^ one starts from (|10|) and rewrites 
the Lagrangian as follows: 

Ltot = L g f + L gh + L(B + a) = 

Lq + L\ + L 2 + L int + Lgf + Lgh (12) 

where Li have the form: 

L 2 {o) = -a v (DlSfj, v - D a b u + igF^a^ = 

= ±aliDl a Dl% u -D™Df-grX>t > ( 13 ) 

Df = d x -6 ca + g f cba B b x = D x 
Lo = -\(F« u (B)f; L x = alD™{B)F^ 

Lint = - l -(D a {B)a v -D u {B)a a y g f abc a\al -\g 2 r hc a\a c J^ a^a{ (14) 

L ghost = -6+(D a (B)D^B + a)) ah 6 b . (15) 
To calculate Green's functions of any hadrons one can use (0) and write 

G h (X;Y) = const J DB uV (B)D^DiljDa^(X)^ in (Y)e Ltot 

= (*-}(X)* in (Y)) Bti , ta . (16) 

Here *&f,i n are initial and final hadron states made of B, a. To calculate the 
integral Da n one can use perturbative expansion in ga^, i.e. neglect in first approx- 
imation L int , containing terms a 3 , a 4 and take into account only quadratic terms 
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L 2 , L g f, L g h. The contribution of L\ was studied in H and it was shown there that 
for the most processes (e.g. for hadron Green's functions) it can be neglected (with 
accuracy of about better than 10%). 

It is convenient to choose background gauge for a M , D^a^ = 0, and take gauge 
transformation in the form 

-> U + (B^ + ~d^U, a M -> U+a^U (17) 

Now one can choose tyf,^ in as gauge- invariant forms built from B^,a^,il),ijj and 
to average over B^ty^^a^ as shown in ([H^). Integrating over DB^ might seem an 
impossible adventure, but we shall see that it can be always written as products of 
Wilson loops with some insertions, which can be treated in two ways: 

i) or in the form of the area law. Then contribution of fields B^ is given and 
fixed by the string tension (r(B). Equating a(B) = a exper one fixes contribution of 
the background fields. 

ii) Using cluster expansion, and expressing result through lowest correlators, e.g. 
(FF), or its scalar parts D(x) and Di(x). In this case background is fixed by these 
functions taken as an input, e.g. from lattice data. 

The rest integrals, over Da^DipD^ are Gaussian and can be simply done. Now 
we turn to the construction of f jin . 

One can choose local or nonlocal forms for ^i n j, the latter are obtained by 
insertion of y) in the local expressions, which are 

*m,/ = $Ji^ for mesons, I\ = 1, 75, 7^, 7^75, o> 

^inj = ^rj/(a, Da)i/j for hybrids, where / is a polynomial, the simplest form is 

^inj — traA 2 a and £r(A 3 a 3 ) for glueballs made of 2 and 3 gluons respectively. 
A 2 ,A 3 are polynomials made of D^B) and ensure the needed tensor structure for 
given quantum numbers. In the simplest case A2 = A3 = 1. 

^in,f = e abc^ ia 1 P b f 2 i34 ! %^K^ ia ^ 2l3 ^ zl for baryons, where a, fi, a are color, flavour 
and Lorentz index respectively. 

It is clear that higher hybrid states for mesons and baryons are obtained by 
adding factors of in ^i n j- 

Insertion of given ^i n j in (0) yields after integration over Da^DipDip the hadron 
Green's function, which can be written respectively for mesons 

G M (X,Y) = (tr(T\ f \x)G q (X,Y)Tt\Y)G ij (X,Y))) B 

- (tr(Ti f \x)G q (X,X))tr(Tt\Y)G q (Y,Y))) B , (18) 
where G q ^ is the quark (antiquark) Green's function, 

G q (X,Y) = @ + m-igB)j c ] Y , (19) 

for hybrids 

G hyb (X,Y) = (tr(T ( i f \x)G q (X,Y)G a (X,Y)Tt ) G q (X,Y))) B , (20) 
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where G g is the gluon Green's function (|TT|). 
For glueballs one has 



G glueball = (tr(A^(X)G g (X,Y)AT l> (Y)G 9 (X,Y))) B (21) 
and for baryons 

G B (X, Y) = (tr(eKW {[ G®(X, Y)eK^)) B (22) 

i=i 

The structure of hadron Green's function is clear: it is a product of quark and 
gluon Green's function averaged over the background field B^. In the next section 
we shall discuss the way, how to make this averaging explicit and express it through 
the Wilson loop. 



3 Feynman-Schwinger representation 

The present stage of development of field theory in general and of QCD in particu- 
lar requires the exploiting of nonperturbative methods in addition to summing up 
perturbative series. This calls for specific methods where dependence on vacuum 
fields can be made simple and explicit. A good example is provided by the so-called 
Feynman-Schwinger Representation (FSR) |9|, |10[ based on the Schwinger proper 
time and Feynman path integral formalism JL1], |12 . 

For some earlier references see |H| and [jT3 
FSR was exploited in [14| and [15| 



More recently some modification of 



The double-logarithm perturbative summation 
of amplitudes is especially convenient for FSR and was extensively studied in [|HJ. 

Meanwhile the first extension of FSR to nonzero temperature field theory was 
done in [|HJ . That was the basis for the systematic study of the role of nonpertur- 
bative (NP) configurations in the temperature phase transitions 
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One of the most important advantages of FSR is that it allows to reduce physical 
amplitudes to the weighted integrals of the averaged Wilson loops. Thus the field 
(both perturbative and NP) enter only through Wilson loops and for the latter 
one can apply cluster expansion theorem E| which allows to sum up a series of 
approximations directly in the exponent. 

Therefore one can avoid Feynman diagram summation to get asymptotics of 
formfactors, as it was done in |L6] . 



The role of FSR in the treatment of NP effects is more crucial. In this case one 
can develop a powerful method of Background Perturbation Theory treating NP 
fields as a background J7], §. 

We start with the scalar particle (e.g. Higgs boson) interacting with the non- 
abelian vector potential via Euclidean Lagrangian 



(23) 
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so that the Green's function of <p can be written using the Fock-Schwinger proper 
time representation as 



poo 

G(x,y) = (m 2 - D 2 X = (x\P J q dse'^^ \y) . (24) 

The next step can be done as in ||, treating s as the ordering parameter. Note the 
difference of the integral (^) from the case of QED treated in [jR], 12, [L3|]: in our 
case is the matrix operator A^x) = y4^(x)T a and it does not commute for different 
x, hence the ordering operator P in (|2"4]). The exact meaning of P becomes clear in 
the final form of the path integral 

roo rx 

G(x,y)= ds(Dz) xy e- K Pexp(ig A^dz^) (25) 

JO Jy 

where K=\ft dr(^) 2 + m 2 s, and 

(^^^n/gl/^i.*.-^., , 26) 



n=l 1 



The last integral in (^) ensures that the path 2^(r),0 < r < s, starts at 
z n{Q) — V/j, an d ends at z^{s) = x^ The form of ( P5|) is the same as in the case of 
QED except for the ordering operator P which gives exact meaning to the integral 
of the noncommuting matrices A^ (zi), A^ 2 (z 2 ) etc. In the case of QCD the form 



(|25"D, ( p6"D was introduced in 0. 

The usefulness of the FSR ( P5|) becomes clear when one considers the physical 
amplitude, e.g. the Green's function of the white state tr(<p + (x)(p(x)) or its nonlocal 
version tr[ip + (x)§(x, y)ip(y)], where $(x,y) to be widely used in what follows is 
parallel transporter along some arbitrary contour C (x, y) 

$(x,y) = Pexp(ig f A^zjdz^). (27) 

One has by standard rules 

G<p(x,y) = (tr(ip + (x)<p(x))tr(ip + (y)ip(y))} A = 

roc roc 

= j Q d Sl J o ds 2 (Dz) xy (Dz') xy e- K - K (W) A +... (28) 

where dots stand for the disconnected part, (G lp (x,x)G tp (y,y))A and we have used 
the fact that the propagator of the charge-conjugated particle ip + is proportional to 

<S> c (x,y) = Pexp(ig f ' Af\z)dz,), A<& = -A™ (29) 
Jy 

and therefore the ordering P must be inverted, &c(x, y) = exp(ig J x A fl (z)dz fl ). 
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Thus all dependence on A^{x) in G v is reduced to the Wilson loop average 

(W) A = (trPcexpig [ A^z)dz^) A . (30) 

Here Pq is the ordering around the closed loop C passing through the points x, y\ 
the loop being made of the path z^ij) and z' (t') to be integrated over in the integral 

Consider now a fermion (quark) Green's function in the Euclidean external glu- 
onic field. Similarly to the QED case considered in |10| one has 

G q (x,y) = (^{x)^{y)) q = (x\(m q + l)) -1 ^) = 
= (x\(m q -D)(m 2 q -D 2 )- 1 \y) = 

/•CO 

(mg-D) ds(Dz) xy e- K ^ a (x,y) (31) 

J 

where ^ a is the same as was introduced in || ^| , 

rX 



x exp(# / dra^F^) (32) 



and a^v = ^(7^7^ — ^vln), while K and (Dz) xy are defined in (|25| ) and fl26|). Note 
that operators Pa, Pf in (|32|) preserve the proper ordering of matrices and a^F^ 
respectively; the explicit examples will be considered below. 

Finally we come to the case of FSR for the valence gluon propagating in the 
background nonabelian field; the detailed study of this case is given in chapter 5; 
here we only quote the result (11) for the gluon Green's function in the background 
Feynman gauge 

G„ u (x,y) = (x\{Dl5^ + 2tgF flu )- 1 \y),D x = d x -tgB\. (33) 

Proceeding in the same way as for quarks, one obtains the FSR for the gluon Green's 
function || 

G^{x,y)= ds(Dz) xy e- Ko <5>^(x,y) (34) 

•J 

where we have defined 

1 r°° ( dz \ 2 f x 
Ko = AJo \~dr) dT] ^^'^ = [ p F p Aexp(ig J A x dz x )x 

x exp(2i<? f dTF ap (z(T))%„ (35) 

J 

Now in the same way as it was done above for scalars in (|2"8|), one can consider a 
Green's function for the physical transition amplitude from a white state of q\, q\ to 
another white state consisting of ^3, q^, i.e the meson Green's function. 

G q q(x, y) = (G g (x, y)TGq(x, y) - G g (x, x)YG q {y, y)) A (36) 

and the first term on the r.h.s. of (|36|) can be reduced to the same form as in ( p8|) but 
with the Wilson loop containing ordered insertions of the operators a^F^ (cf.(|32|)). 
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4 FSR at nonzero temperature 



This part of lecture is based on papersO, [18|, where FSR was first introduced for 



T > 0. We discuss first basic formalism at T > and then gluon and quark Green's 
functions. 

4.1 Basic equations 

We start with standard formulas of the background field formalism f?|, |] generalized 
to the case of nonzero temperature. We assume that the gluonic field can be 
split into the background field and the quantum field a M 

A„ = B fl + a„, (37) 

both satisfying periodic boundary conditions 

B^(z A , Zi) = B^(z A + n/3, Zi), 

a^(z 4 , = a^(z 4 + nP, z { ), (38) 
where n is an integer and j3 = 1/T. The partition function can be written as 

Z(V,T,n = 0) = (Z(B)) B , 
Z(B) = N J D0exp (- £ 1 dr J d 3 xL tot (x, r)^J , (39) 

where denotes all set of fields a M , N is normalization constant, and the sign 

( )b means some averaging over (nonperturbative) background fields B^, the exact 
form of this averaging is not needed for our purposes. Here L tot is the same as in 
(12)-(15). 

Integration over ghost and gluon degrees of freedom in (|3§|) yields 
Z{B) = iV'(det W(B)); e f[det(-D,(B)D,(B + a))] a= s x 

o J 

x ll + ±?f(a= £)}exp(- \JW-*J) • (40) 

I 1=1 11 \ dJ J) \ 2 / J»=D^(B)F^(B) 

One can consider strong background fields, so that gB^ is large (as compared to 
Aq CD ), while a s = g 2 '/An in that strong background is small at all distances. More- 
over it was shown that a s is frozen at large distances |7], |8], |2Efl . 

In this case ( ^0|) is a perturbative sum in powers of g n , arising from expansion 
in {ga^y. 

In what follows we shall discuss the Feynman graphs for the free energy F(T), 
connected to Z(B) via 

F(T) = -T\n(Z(B)) B . (41) 

As will be seen, the lowest order graphs already contain a nontrivial dynamical 
mechanism for the deconfinement transition, and those will be considered in the 
next section. 
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4.2 The lowest order gluon contribution 

To the lowest order in ga^ (and keeping all dependence on gB^ explicitly) one has 

Z = e - F ^' T = N'(exp(-F (B)/T)) B , (42) 
where using fllPf) F (B) can be written as 

-F (B) = ^\ndetW -\jidet(-D 2 (B)) = 



\(t)-e- tw + /°°C(t)-e 4D2(B) l . (43) 
t Jot] 



In ( |4"3"| ) implies summing over all variables (Lorentz and color indices and 

d M 2a t s 
ds V(s) 



coordinates), ((t) = lim 4- ^/f | s =o is regularizing factor, one can use also the Paul! 



Villars form for C{t)- 

Graphically, the first term on the r.h.s. of fl4"3"| ) and is a gluon loop in the 
background field, the second term is a ghost loop. 

Let us turn now to the averaging procedure in (42|). With the notation tp = 



—Fq(B)/T, one can exploit in (f42|) the cluster expansion [^2 

(exp^) B = exp (£<<^>>^J = exp{<^) B + -[<#)% - tf) B ] + O(^)}. (44) 

To get a closer look at (<p)b we first should discuss thermal propagators of gluon 
and ghost in the background field. We start with the thermal ghost propagator and 
write the FSR for it [g 

(-D 2 )- 1 = (x\ I dte w2 ^\y} = j o dt(Dzr xy e' K Hx,y). (45) 

Here $ is the parallel transporter in the adjoint representation along the trajectory 
of the ghost: 

$(z,y) = Pexp(^ y Bptydzp), (46) 

also 

K= lh drZ » Z » = — 
and {Dz)™ y is a path integration with boundary conditions imbedded (this is marked 
by the subscript (xy)) and with all possible windings in Euclidean temporal direction 
(marked by the superscript w). 
One can write it explicitly as 



N d 4 ((m) d 4 p [ / N N 



(47) 
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Here, ((k) = z(k) — z(k — 1), Ne = t. 

One check that in the free case, B^ = (|45|) reduces to well-known form of the 
free propagator 



,2\-l 

I XI) 



dtexp 



i 



45 



d 4 p 



with 



x exp ]ip (J2 CM - - V) ~ nPSfjAj] 

poo r 

E / exp — p 2 t — — z) — ip^nP 

n J° 

d((m) 



dt 



(2tt)< 



(48) 



d((m) 



(47T£) 



Using the Poisson summation formula 



IT- E exp(ip 4 nf3) = E S(p^-2nk) 

A7T n=0.±l,±2... fe=0.±l,... 



one finally gets the standard form 



-0 



2X-1 

l.ilj 



E 

fc=0,±l, 



Td 3 p exp[— ipi(x — y\ — i27rkT(x 4: — y 4 )] 



(27T) 



p? + (2ttA;T) 2 



(49) 



(50) 



Note that, as expected, the propagator (|45|),(|5"0"D corresponds to a sum of ghost 
paths with all possible windings around the torus. The momentum integration in 
( f4~7|) asserts that the sum of all infinitesimal "walks" C( m ) should be equal to the 
distance (x — y) modulo N windings in the compactified fourth coordinate. 

For the gluon propagator in the background gauge we obtain similarly to fl4"5]) 



dt(Dz)™ y e- K $ F (x,y), 



where 



$ F (x,y) = PfP^xp [2ig F(z(r))dr) exp [ig \ B^dz^ 



(51) 



(52) 



and the operators PfP are used to order insertions of F on the trajectory of the 
gluon. 

Now we come back to the first term in fl44D , (<p) b, which can be represented with 
the help of (gSp and (|1]) as 



B 



((t)d 4 x(Dz)~ x e 



-tr(<& F (x, x))b — (tr$(x,x)) 



B 



(53) 



where the sign tr implies summation over Lorentz and color indices. 

One can easily show (see || for details) that (|53| ) yields for £> M = the usual 
result for the free gluon gas: 



4^2 



F (B = 0) = -T<p(B = 0) = -(iV c 2 - 1)V 3 



~l5 



(54) 
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4.3 The lowest order quark contribution 

is done trivially and leads to the following 



Integrating over quark fields in 
additional factor in fl40|) 



det(m + D(B + a)) = - det(m 2 - D 2 (B + a)), 



(55) 



where we have used the symmetry property of eigenvalues of D. In the lowest 
approximation, we omit a M in ( |55| ) and write the free energy contribution Fq(B) 
similarly to the gluon contribution (f42"D,(|43D as 



±F 9 (5) = --lndet(m 2 - D 2 {B)) = --Sp / C(t)-e~ tm ' +tL ' (B) , 



dt 



D'- 



(56) 



{D^f = Dl^+gF^a^^D 2 + gaF; ^ = -(7^-7^). (57) 



T 

where the sign Sp has the same meaning as in (fj3|) and 



Our aim now is to exploit the FSR to represent ( [56D in a form of the path integral, as 
it was done for gluons in ([25]) . The equivalent form for quarks must implement the 
antisymmetric boundary conditions pertinent to fermions. It is easy to understand 
that the correct form for quark is 



where 



^(B) = -~tr jTc(*) jd±xjDzr xx e- K - tm2 W a (C n ), 
W a (C n ) = P F P A exp (ig A^dz^j exp g (aF) , 



(5? 



(Dz] 



xy 



1 (4^) 2 n= o,±l,±2, 



(2tt)< 



exp 



■ / N 

ip E c( 

Vm=l 



m) — (x — y) — n/3<5 



(59) 



One can easily check that in the case B^ = one is recovering the well known 
expression for the free quark gas 



7tt 2 



F 9 (free quark) = -—N C V 3 T ■ n f 
loll 



(60) 



where nf is the number of flavors. The derivation of ( |5"0"D starting from the path- 
integral form (|58D is done similarly to the gluon case given in the Appendix of [§. 
The loop C n in (|5*%|) corresponds to n windings in the fourth direction. Above the 



deconfmement transition temperature T c one can visualize in (|58|) the appearance 
of the factor 



Q = exp 



ig \ B i {z)dz 4: 
Jo 



(61) 
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(Note the absence of the dash sign in (pi] ) as compared to(^) implying that in 
( |B~TD B 4 is taken in fundamental representation). For the constant field B4 and 
B; — 0, i — 1, 2, 3, one obtains 



n=l 



\n+l 



n 1 



(62) 



This result coincides with the obtained in the literature |27 



5 Wilson loop and Field Correlators (FC) 

Hadron Green's functions in the form of path integrals due to FSR are expressed 
through Wilson loops; moreover all gluonic fields enter via these loops, and in case 
of quarks and gluons additional spin factors appear, namely 

W F (C) = trP A expig f ■ P F exp f g ]T F(z(r))dr (63) 

J Co J 

where J2 F 1 — o^F^ for quarks and 2iF flu for gluons, and the dash sign marks the 
adjoint representation for gluons in contrast to the fundamental one for quarks. 



As it was shown previously in P, B3, [28 , the spin factors can be treated per- 



turbatively in the effective Hamiltonian and this topic will be discussed in the last 
lecture. Here we shall consider the simplest form of Wilson loop in the arbitrary 
representation D of the SU(N), 

W(C)=tr D P A exptg [ A^, A fl = A a T a . (64) 

J Co 

The nonabelian Stokes theorem |2JJ allows to express W(Cq) in terms of field 
strength F^, instead of the vector potential A^, this later will enable us to express 
the vacuum average (W(C )) through field correlators. 

The easiest way to derive the Stokes theorem is via the generalized coordinate 
gauge |}0, ^H| , namely one introduces the following representation for A^ 

A,(x) = £ 8z f x ' x) dz u F vp {z{s, x)) (65) 

Here we have introduced an arbitrary contour C, which is defined by a coordinate 
z(s, x), starting at z = x,z(l,x) = x, while for s = 0,z(0,x) = xq;xo may also lie 
at infinity. 

The only condition which one should impose on the contour C(x) is the "auto- 
morfic condition" [pi], i.e. 

zJs, x) = zJs", z(s', x)) (66) 
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which means that choosing on the contour C(x) some point x' = z(s',x) to identify 
it with the initial point of the new contour C(x'), one gets C(x') lying entirely on 
the old contour C(x). 

Inserting now ((65|) in (|64]) one obtains the nonabelian Stokes theorem 



W(C )=tr D P s expig f do pv (z)F pv {z) (67) 

JS(Cn) 



'S(C ) 

Note that the ordering P$ in (|67D is fully determined by the contour C(x) in 
and ordering of points x along the contour Cq. 



One can show |3T| that parallel transporters ^i(x,y) = P expig J L A^dz^ in 
the generalized coordinate gauge are equal to unity whenever contour L(x,y) lies 
on the contout C(x). Therefore one can render the Stokes theorem ( |67D to the 



gauge- covariant form 

W(C Q ) = tr D P s expig do^^F^z^Xo) (68) 

J s 

where 

F^z, x ) = $(x , z)F flu {z)<5>{z, x ) (69) 

The crucial step now is to use the cluster expansion theorem [22 and to retain 
only the lowest term, namely 

(W(C )) = (tr D P s expig J^da pi/ (z)F pu (z, x )) = 

= tr D exp j-y J^J^da^ u (z)da pu (u){{F^(z,xo)F pa (u,xo))) + ...j . (70) 

Here dots imply higher order terms, e.g. ((FFF)), etc. Moreover double brackets 
mean cumulants, e.g. ((AB)) = (AB) — (A)(B). The angular brackets imply 
averaging over all gluonic fields with the standard QCD action, as in (TO). 

Next step is the standard representation of the bilocal field correlator in (|70|). 
made in |32| . It goes as follows. First one realizes, that when xq is somewhwre in the 



middle of surface S, then distances \z — Xq\, \u — x \ are of the order of R - "radius" 
of S, while average distance \z — u\ will be of the order of the vacuum correlation 
length T g . 

Hence in the situation when i? 3> T g (and this is typical situation for hadrons, 
since as we shall see R ~ size of hadron ~ 1 fm, T g ~ 0.2 fm) one can approximate 

((F^(z, x )F pa (u, x )» = (F^(z)$(z, u)F p(J {u)<b{u, z)) + O ( ) (71) 



The lowest correlator ( |7T| ) with parallel transporter connecting points x and y 
along the straight line has a general representation in terms of two scalar form-factors 
D and D 1 W% 

1 

—tr(F pu (x)$(x,y)F Xa (y)$(y,x)) = 
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= (<V<W - S fia S v x)D(x -y) + -d^hxSua - h a 5 vX )+ 

+ perm]D 1 (x - y)}, h il = x ii -y tl . (72) 

It is important here that by definition the term with D\ is a full derivative. In 
what follows both functions D and Di will have a basic role in all nonperturbative 
dynamics in QCD, both at zero and nonzero temperature. To understand why this 
role is so important one must compare contribution of the bilocal FC ( |72|) with all 
higher correlators, neglected in (|70|). 

That comparison can be done at present only indirectly since the direct deter- 



mination of higher FC on the lattice is difficult and still not accurate enough |33| . 

Happily there was another type of lattice measurement recently ]33| which has 
allowed to estimate the contribution of higher FC with unprecendented accuracy 
-better then one percent! 



Since the corresponding work [go] is important for understanding the structure 
of the QCD vacuum, we shall discuss it is some detail. 

Consider static sources in representation D of SU (3) and define static potential 
between sources through the Wilson loop as follows 

V D {R) = - lim \HW{C)) (73) 
where the Wilson loop in the 34 -plane has the form RxT and cluster representation 

oo „ 

(W(C)) = Tr D exp £ / (ig) n ((F(l)...F(n)))da(l)...da(n). (74) 

n=r] 

Here we have omitted Lorentz indices and denoted 

F(n)da{n) = F{x {n \x ) = $(x , x {n) )El(x {n) )T a <S>(x {n \ x )da u {x) (75) 

The SU(3) representations D = 3,8, 6, 15a, 10, 27, 24, 15s are characterized by 
3 2 — 1 = 8 hermitian generators T a which satisfy the commutation relations 
[T a ,T b ] = if abc T c . One of the main characteristics of the representation is 
an eigenvalue of quadratic Casimir operator C D , which is defined according to 

»(2) 



C D = T a T a = Cd ■ 1- Following the notations from |34| we introduce the Casimir 



ratio dr, = Cd/Cf, where the fundamental Casimir CV = (AT? — 1)/2N C equals to 
4/3 for SU(3). The invariant trace is given by Tr^l = 1. 

Since a simple algebra of the rank k has exactly k primitive Casimir-Racah 



operators |36[ of order m 1 , ..,m^, it is possible to express those of higher order in 
terms of the primitive ones. In the case of SU(3) the primitive Casimir operators 
are given by 

d§ = 5 ab T a T b ; Cg> = d abc T a T b T c (76) 



while the higher rank Casimir operators are defined as follows 

j(r) 

\h-ir) 



c£ = d ( £; ) T i \.T i r (77) 
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where the totally symmetric tensor dffi ^ on the SU(N C ) is expressed in terms of 
5ik and d^ k (see, for example, |37|). 

The potential ( |73|) with the definition ( [73] ) admits the following decomposition 

V D (R) = d D V {2 \R) + d 2 D V^\R) + (78) 

where the part denoted by dots contains terms, proportional to the higher powers 
of the quadratic Casimir as well as to higher Casimirs. 

The fundamental static potential contains perturbative Coulomb part, confining 
linear and constant terms 

V D {R)=a D R-v D - e -§ (79) 



The Coulomb part is now known up to two loops |38] and is proportional to Cd- 



The "Casimir scaling hypothesis" [159] declares, that the confinement potential is 
also proportional to the first power of the quadratic Casimir Cd, i-e. all terms in 
the r.h.s. of (|78|) are much smaller than the first one. In particular, for the string 
tensions one should get (Jd/cf = do- 

This scaling law is in perfect agreement with the results found in [|4|, . Earlier 
lattice calculations of static potential between sources in higher representations 
are in general agreement with [|34] . 



To see, why this result (to be more precise - why the impressive accuracy of the 
" Casimir scaling" behaviour) is nontrivial, let us examine the colour structure of a 
few lowest averages in the expansion (|71|). 

The first nontrivial Gaussian cumulant in ([F4]) is expressed through Cd and 
representation-independent averages as 

Tr D (F(l)F(2)> = (F«(l)F a (2)} = ^ (F«(l)F a (2)), (80) 

C c 

so Gaussian approximation satisfies " Casimir scaling law" exactly. It is worth being 
mentioned, that this fact does not depend on the actual profile of the potential. It 
could happen, that the linear potential observed in |3JJ] is just some kind of inter- 
mediate distance characteristics and changes the profile at larger R (as it actually 
should happen in the quenched case for the representation of zero triality due to the 
screening of the static sources by dynamical gluons from the vacuum, or, in other 
words, due to gluelumps formation). The coordinate dependence of the potential 
is not directly related to the Casimir scaling, and can be analized at the distances 
which are small enough to be affected by the screening effects. 

Having made these general statements, let us come back to our analysis of the 
contributions to the potential from different field correlators. We turn to the quartic 
correlator and write below several possible colour structures for it. We introduce 
the following abbreviation 

/pM) — (F a (l)T a F b (2)T b F c (3)T c F d (4:)T d ) — j ia T b T c T d (F^) abcd 
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where the Lorentz indices and coordinate dependence are omitted for simplicity of 
notation. One then gets, with some work in the last case the following possible 
structures 

(FWf bcd ~ 5 ab 5 cd (FW) ~& D .l 

(FW) abcd ~ 5 ac 5 bd (FW) ~ (CI - \N C C D ) • 1 

(F^ ~ fadefd* (F [4] ) ~ -\N 2 C C D ■ 1 

(F^ ~ facef** (F^) = 

(FW) M ~ f apm f bpn f dem f cen (F®) ~5(CS+ |C D ) ■ 1 

(3) (2) 

The operator Q/ enters together with C}/ at higher orders. Notice, that the 
terms, proportional to the square of Cd appear in both the 55 parts (the first and the 
second strings) and higher order interaction parts (the last string). Mnemonically 
the Cd - proportional components arise from the diagrams where the noncompen- 
sated colour flows inside while the C D - components describe the interaction of two 
white objects. It is seen, that the Casimir scaling does not mean "quasifree gluons", 
instead it means roughly speaking "quasifree white multipoles" (see discussion at 
the end of the paper). 

Let us analyse the data from |M| quantitatively. We have already mentioned, 



that the Coulomb potential between static sources is proportional to up to the 
second loop (and possibly to all orders, this point calls for further study) and hence 
we expect contributions proportional to C D ~ d D to the constant and linear terms, 
i.e. we rewrite ([79]) as follows 

V D (R) = d D vV\R) + dUvff + ct { d } R). (81) 

and all higher contributions are omitted. Here v^^a^ measure the (^-contribu- 
tion of the cumulants higher than Gaussian to the constant term and string tension 
respectively. The results of the fitting of the data from [34] with (]81f) for some 
representations are shown in the Table 1. 

All numbers in the Table 1 are dimensionless and given in lattice units. The 
author of p3] used anisotropic lattice with the spatial unit a" 1 = 2AGeV. 

Several comments are in order. First of all it is seen that the Casimir scaling 
behaviour holds with very good accuracy, better than 1% in all cases in the table 
1 with the reasonable x 2 . It should be stressed, that any possible systematic errors 
which could be present in the procedure used in [[34] must either obey the Casimir 
scaling too or be very tiny, otherwise it would be unnatural to have the matching 
with such high precision.]] Nevertheless, the terms violating the scaling are also 
clearly seen. While the value of the constant term v$ is found to be compatible 
with zero within the error bars, it is not the case for . We have not found any 



1 Since we are mostly interested in the relative quantities, their actual magnitude in the physical 
units is of no prime importance for us. This is another reason why we do not discuss possible 
systematical errors of [B4[ and finite volume effects. 
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Table 1: The Casimir-scaling and Casimir-violating string tensions and shifts. 
Based on the data from G.Bali, |hep-lat/990802Tl . [34] 



D 


ag } • 10 4 


Aag } • 10 4 


v$ ■ 10 4 


Aug 5 • 10 4 




X 2 /N 


8 


-3.486 


1.2 


-2.513 


2.8 


0.004 


19.22 / 43 


6 


-6.428 


1.2 


0.950 


2.6 


0.007 


25.76 / 42 


15a 


-5.244 


0.55 


-0.5611 


1.1 


0.003 


39.06 / 42 


10 


-4.925 


0.50 


0.2489 


1.0 


0.003 


22.05 / 41 



sharp dependence of a D on the representation D, which confirms the validity of 
the expansion (|8~I| ) and shows, that the omitted higher terms do not have significant 
effect in this case. Notice the negative sign of the string tension correction. In 
euclidean metric it trivially follows from the fact, that the fourth order contribution 
is proportional to (ig) 4 > while the Gaussian term is multiplied by (ig) 2 < for 
real g. 

The analysis given above allows to draw important conclusions about the QCD 
vacuum: 

1) The contribution of higher FC to the static potential is less than 1% at all 
measured distances 0.1/m < R < 1.3/m. (Another possible scenario is that all 
higher FC are fine-tuned to cancel exactly all terms proportional to C£>, n > 1; but 
this fine-tuning seems rather artificial). 

2) The lowest correlators D and D\ should describe dynamics of the Wilson loop 
(i.e. the dynamics of nonperturbative QCD with accuracy better than 1%. This fact 
can be checked in calculations of physical amplitudes, e.g. we shall give exmples of 
our predictions in the last lecture. 

3) All modes of the QCD vacuum which are based on classical solutions are 
strongly limited by the observed Casimir scaling, since a classical lump (e.g. an 
instanton or a dyon) generate a big Quartic (and higher) FC. 

E.g. the density of instanton gas (in the SU(2) case) in the QCD vacuum has an 
upper bound of y — 0.2/m -4 (and even smaller if one considers smaller distances). 
This should be compared with the normal density of instanton gas of lfm~ 4 pro- 
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viding the standard gluonic condensate. 

The picture of the QCD vacuum which emerges from the Casimir scaling is the 
picture of stochastic ensemble of small white dipoles, where only any two adjacent 
points are correlated, while all the rest are statistically independent. This picture 
is ideal for description in terms of lowest field correlators, D(x) and D\{x) at all 
temperatures, and we shall indeed show how easily one can describe dynamics of 
different precesses in QCD. 

4) As a by-product the results of [34,35] give a strong upper bound on the so- 
called Luescher-Weissz term in the potential f40|, Vlw = — jf^, since the latter 
violates Casimir scaling. This means that the QCD string is either not the Nambu- 



Goto string, used in |4(|, or else one should use properly the finite lattice size for 
vibrating strings - this is an interesting problem for the future. 



6 Field correlators at zero and nonzero tempera- 
ture 

Since D(x) and D\(x) are basic elements of the nonperturbative QCD dynamics, it 
is important to find them as functions of distance x and use later on as an input. 

The only standard way to get information on D, Di before was the lattice sim- 
ulation method, and lattice measurements of D,Di have been done both at T = 



IJ at T > |2j. 

At the end of these lectures we shall describe other methods of calculating field 
correlators and in particular D and Di, but now we shall describe in detail results 
of the Pisa group. 

For T = these data are presented in Fig.l, and can be parametrized by the 
form 

D(x) = Ae» x + (82) 



x 4 



Up Mi 1 

D 1 {x) = Be» x + — -j- (83) 



;r 4 



The last terms in (^2|), ( |S3l) proportional to x~ A are clearly of perturbative origin, 
moderated at large distances by screening factors. The most important information 
is contained in the first terms in ([82]), (|83|) proportional to exponents. It is striking 
that numerically \i — // ~ 1 GeV for the unquenched case. 



At the end of the lectures we shall explain the appearance of exponents in (82 ) . 
( |8~3"D and the mass values around 1 GeV - this will be connected to the masses of 
the so-called gluelumps - bound states of gluons with a static adjoint source. These 
masses have been calculated recently both on the lattice |^3[ and analytically 



Physically appearance of large masses - of the order of 1 GeV leads to very 
important consequences: 
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i) First, the gluon correlation length T g = 1/// « 0.2 /m is small as compared to 
the generic hadronic size and implies that stochastic picture of small independent 
white dipoles which was discussed in previous section. 

ii) As a consequence of small T g one can develop the potential picture of hadrons 
(see |23], [28| for more discussion), with quarks and valence gluons interacting via a 
string instantaneously without retardation, while nonpotential corrections are pro- 
portional to some powers of (T g /R h ), where Rh is the typical hadron size. 

iii) Finally, smallness of T g puts serious limitation on the use of some standard 
methods in QCD. E.g. the application of QCD sum rules strictly speaking implies 
large values of T g , namely T g ^> Rh, otherwise the notion of gluon condensate and 
its derivatives loses its meaning. 

The same is true for the Voloshin-Leutwyler method, based on the constant 
condensate contribution to the heavy quarkonia dynamics. One can estimate that 
the inequality Rqq < T g holds true for toponium and hardly for bottomonium. 

Knowing D and D\ one can calculate all the dynamics of heavy quarkonia (and 
also of other hadrons - in the next sections). 

In particular the string tension is defined by the nonperturbative part of D(x) 



a 



^Jd 2 xD(x) (84) 



Note that D\ does not contribute, being the full derivative, and the perturbative 
part of D(x) is cancelled by the corresponding perturbative parts of triple and higher 
correlators ||45|| . The whole static potential, as expressed through D and D\ has the 
form 



rR roo rR roo 

V{R) = 2R dX dvD(X,v) + XdX du[-2D(X,v) + D 1 (X,v)] (85) 
Jo Jo Jo Jo 

Here D(X, v) = D(\JX 2 + u 2 ) because of Euclidean 0(4) invariance. Note that the 
first term on the r.h.s. of ( jS5|) yields at large R the term aR, while the second term 
adds there a constant. The perturbative contribution to V(R) can be included in 
D\, or considered separately, in both cases one obtains the form 

V{R) = aR - v - 4/3^4^ (86) 

R 



observed many times on the lattice and recently in [[34] 



We shall comment in what follows on the NP dependence of a s (R), but now 
we turn to spin-dependent corrections to the potential, calculated in |2^, |46|, i.e. 
spin-orbit potentials Vx, V2, tensor force V3 and spin-spin potential V4 (we exploit the 
standard notations of Eichten and Feinberg j|7||). All these potentials are expressed 



through D and D\ only and in this way the dynamics of heavy quarkonia is fully 
determined. 

The corresponding calculations have been done repeatedly and successfully dur- 
ing last 10 years in HH-Hi. In the Tables 2 and 3, given are the known levels of 
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charmonium and bottomonium calculated using the lattice data on P, Pi from in a 
standard way. One can see a very good description of all levels, which far exceeds 
in quality the newly developed NRQCD methods. 

In the Table 2 listed are bottomonium level differences, which have more physical 
significance than the absolute values of levels, since the heavy quark mass is not 
known and used as an input parameter. Calculations in [51] have been done both 
in relativistic and nonrelativistic kinematics, and one can estimate therefore the 
accuracy of 1/M expansion usually done in the framework of NRQCD. One can 
see that corrections can be very significant, especially for P-wave levels, where the 
matrix element of the 1/r 3 operator is increased up to 20% in relativistic calculations 



0- 

Similar calculations have been done for charmonium in |)U] , again with very good 
results. An interesting comparison of these and similar results is contained in (52J. 
where the Voloshin-Leutwyler method was used, and shown to be applicable (if any) 
only to ground state of bottomonium. 



Table 2: Bottomonium level differences (MeV) for the Schrodinger and Salpeter 
equations from [jSlj] 



Masses differences 


Set I, a 


eff = 0.3545 


Set II, a e ff = 0.36 


Exp.val 






m = 4.737 GeV 


m = 


4.81GeV 


(MeV) 






a = 0.20 GeV 2 


a = 0.18 GeV 2 






Rel. 


Nonrel. 


Rel. 


Nonrel. 




M(2S) 


- M (15) 


554.34 


551.97 


556.55 


550.03 


562.9±0.5 


M(3S) 


- M (25) 


350.43 


354.78 


335.62 


338.49 


332.0±0.8 


M{4£) - 


■M(3S)( a ) 


285.93 


291.83 


270.63 


275.30 


224.7± 4.0 


M(1P) 


- M(1S) 


458.04 


439.66 


473.49 


450.15 


439.8±0.9 


M(2P) 


-M(1P) 


359.67 


366.75 


342.55 


348.70 


359.8±1.2 


M{2S) 


- M (LP) 


96.31 


112.31 


83.07 


99.88 


123.1±1.0 


M(3S) 


- M(2P) 


87.06 


100.34 


72.82 


89.67 


95.3±1.0 



( a ) The 4S level lies above the bB threshold and the corresponding threshold 
correction was not taken into account, which explains the discrepancy. 
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Table 3: The spin-averaged mass level differences in charmonium a ) (in 
MeV) 





Experimental 
values 6 ) 
(MeV) 


SetC 
M = 1.48 GeV, 
a = 0.18 GeV 2 , 
a s {no) = 0.365, 
^ = 0.909 GeV, 
a v = 0.42 


Set D 
m = 1.40 GeV, 
o = 0.183 GeV 2 , 
a fl (/x ) = 0.312, 
fi = 1.357 GeV, 
a v = 0.39 


M(2S) - M(1S) 
M(1P) -M(1S) 
M(2P) - M (IP) 


595.39 ± 1.91 
457.92 ± 1.0 
137.47±1.77 


591.30 
460.68 
130.62 


586.09 
445.93 
140.15 



a ) all matrix elements defining spin effects were calculated for spinless Salpeter 
equation. 

6 ) M(nL) means the spin-averaged mass. 

7 The Field Correlator Method for light quarks 

Let us start with the meson Green's function and use for that FSR for both quark 
and antiquark. 



G M (x,y) = (trT^(m-D) ds dse~ K ~ K (Dz) xy (Dz) xy T { - i \m — D)W F ) (87) 

J « 

Here the bar sign refers to the antiquark, and 

W F = P A P F exp(z 5 j dz^) exp(g j* dra^F^) exp(-g jT a^F^df) (88) 



In (|87D enter integrations over proper times s, s and r, f, which also play the role of 
ordering parameter along the trajectory, = z^t), z^ = z^(f). 

It is convenient to go over to the actual time t = z± of the quark (or antiquark), 
defining the new quantity //(£), which will play very important role in what follows 



2fi(t) = ^-, t = z 4 (r) 
dr 
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For each quark (or antiquark and gluon) one can rewrite the path integral ( |54"D , (|8~TD 
as follows 

jf ds(D A z) xy ... = const J Du{t)(D*z) xy ... (90) 
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where (D 3 z) xy has the same form as in fl2"E|) but with all 4- vectors replaced by 3- 
vectors, and the path integral Dfi(t) is supplied with the proper integration measure, 
which is derived from the free motion Lagrangian. 

In general can be a strongly oscillating function of t due to the Zitterbewe- 
gung. In what follows we shall use the steepest descent method for the evaluation 
of the integral over D/j,(t), with the extremal fio(t) playing the role of effective or 
constituent quark mass. We shall see that in all cases, where spin terms can be 
considered as a small perturbation, i.e. for majority of mesons, fiQ is positive and 
rather large even for vanishing quark current masses m, m, and the role of Zitterbe- 
wegung is small (less than 10% from the comparison to the light-cone Hamiltonian 
eigenvalues, see |53|, |54| for details). 

Now the kinetic terms can be rewritten using (^) as 

K+R - + %wr + 11 + 2I) + ^r |(i - (t))2 + 11} (91) 

where T = x 4 — 2/4. In the spin-dependent factors the corresponding changes are 



dt 

--a*-^ j Q 2 ^ 



dra^F^ = / —— a ^F^(z{t)). (92) 



In what follows in this section we shall systematically do perturbation expansion of 
the spin terms, which contribute to the total mass corrections of the order of 10- 
15% for lowest mesons and much smaller for high excited states. This perturbative 
approach fails however for pions (and kaons) where chiral degrees of freedom should 



be taken into account. In this case other equation should be considered |p5 
Therefore as the starting approximation we shall use the Green's functions of mesons 
made of spinless quarks, which amounts to neglect in (p7D,(p8|) terms (m — D),(rh — 



D) and o^F^. As a result one has 



G$(x,y) = const J D^t)Djx{t){D z z) xy {D z z) xy e- K - R \W). (93) 

Our next approximation is the neglect of perturbative exchanges in (W), (they 
will be restored in the final expression for Hamiltonian) which yields for large Wilson 
loops, .R, T > T g , 

(W) B = const exp(-aS min ) (94) 
where S m i n is the minimal area inside the given trajectories of quark and antiquark, 

Smin = J dt dfi^detg, g ab = d a w lx d h w' x , a, b = t, (3. (95) 

The Nambu-Goto form of S m i n cannot be quantized due to the square root and 
one has to use the auxiliary field approach [^7[ with functions v(/3,t) and n((3,t) to 
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get rid of the square root, as it is usual in string theories. As the result the total 
Euclidean action becomes |58| 



A = K + K + aS mm = dt d[3{-(— + —) + ^^R 2 + 

Jo Jo 2 n(t) p,(t) 2 



f + ^[^ 2 + (-)V 2 - 2r/Or) + ?? 2 r 2 ]}. (96) 



Here fi + = li + fx, fx = J ^=, R4 = , r , = z, - z». Note, that integrations 

over 11, v and rj effectively amount to the replacement by their extremum values |58| . 

Performing Gaussian integrations over and 77 one arrives in the standard way 
at the Hamiltonian (we take m = fa for simplicity) 

_p 2 + m 2 L 2 /r 2 

~ ~WV + m + m+ zfoW -±)m/w + 

^f* + [ 1U ^, (97) 



/o 1/09) 

where p 2 = (pr) 2 /r 2 , and L is the angular momentum, L = (r x p). 

The physical meaning of the terms n(t) and z/(/5) can be understood when one 
finds their extremal values. E.g. when a = and L = 0, one finds from 



^ = 2 V /p 2 + m 2 , / i = V /p 2 + m 2 (98) 

so that liq is the energy of the quark. Similarly in the limiting case L —>■ 00 the 
extremum over v{[3) yields: 



or 



"o(/?)= t f= f= , H 2 = 2na^/L(L + 1) (99) 



so that Uq is the energy density along the string with the (3 playing the role of the 
coordinate along the string. 

Let us start with the L = case. Taking extremum in v{(3) one has 

2 2 

Li[t) 

Here /x(t) is to be found also from the extremum amd is therefore an operator in 
Hamiltonian formalism. Taking extremum in fi(t) one obtains 



H 2 = 2^/p 2 r + m 2 + ar. (101) 

The Hamiltonian ( |101|) is what one traditionally exploits in the relativistic quark 
model (RQM) |5£| (apart from color Coulomb and spin-dependent terms to be dis- 
cussed below). The RQM was an essential step in our understanding of hadronic 
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spectra. At the same time the usual input of RQM contains too many parameters 
and the model was introduced rather ad hoc. Another deficiency of this model at 
this point is twofold: 

i) one usually takes m in (|9S|) , ( |101 ) to be constituent quark mass of the order of 



100-200 MeV, which is introduced as an input. Instead we have in (|97| ) ,(|9^ ) ,( |101| ) 
the current quark mass renormalized at the reasonable scale of 1 GeV. Hence it is 
almost zero for light quarks; 

ii) the form QlOlj) is used in RQM for any L, and to this end one writes in ( |101|) 

vl p 2 - 

However the Regge slope for both ( |100| ) and ( |101| ) is l/8cr instead of the string 
slope which occurs for the total Hamiltonian (p7|) , since the RQM Hamiltonian 
( |101| ) does not take into account string rotation. 

Still for L = the form ( |101| ) is a good starting approximation, and it is reward- 
ing that our systematic approach makes here contact with RQM. Sometimes it is 
convenient to use instead of (|101| ) a more tractable form ( |100|) where = fi and 
/io is to be found from the extremum of the eigenvalue of Hamiltonian H\\ 

H x tl) = M x il), M 1 =Mi(ji ), ^ = 0. (102) 

dfj, 

Equation (|100|) , QlOlj) can be reduced to the Airy equation with the eigenvalues 
written as (for m = 0) 

M«=4/i (n), /xo(n) = ^(^f 1 ' (103) 

where a(n) is the corresponding zero of the Airy function. The first comprehensive 
study of the equation (|102|) with fixed fi(t) = m was done in where also a(n) 
are quoted. 



It is easy to apply the Hamiltonian ([37]) to the case of light, []23], |6T| heavy [| 



[ plR and heavy-light mesons |24], |62fl. We refer here to where these applications 
and also hybrid calculations have been successfully done. Again, the most important 
point of the described formalism is the natural occurence of the constituent mass 
of quark and gluon, which is calculated through the string tension (and a s ) in the 
unique way, as shown in Eq.( |102| ). 



Since this point is very important from both phenomenological and physical 
point of view, we shall discuss here in more detail the case of baryon magnetic 
moment (defined through quark constituent mass) and spectrum of glueballs (where 
constituent mass of gluons enter in a decisive way). 

For baryons the Hamiltonian is derived from the Green's function of a baryon in 
the same way as it was done above for mesons. 

We refer the reader to fl53]-||65|| for details of derivation. The result is 
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Here m k are current quark masses defined at the scale of 1 GeV, fi k - the corre- 
sponding constituent masses (energies) to be calculated in what follows, r( fc ) is the 
vector from position of the k-th quark to the string junction. The coordinates r\ 
are Jacobi combinations defined through in a standard way. 

The eigenvalues of (HMD Mjji) can be found in a simple way with good accuracy 
using the the hyper-spherical explansion method (menthod of the K - harmonics) 
If)!! . Assuming equal masses M k = m and hence equal constituent masses = fi, 
one obtains equation for the definition of /j, 



dM n (fj) 



/,From that we find fi = fi u for baryon consisting of 3 light quarks 

o / o \n 3 /4 



. (105) 




3 • 5V3 I 375 



0.957^a. (106) 



Choosing a =0.15 GeV 2 (as it accepted for baryons phenomenologically |67j and 
recently found microscopically p5[ ) one has \x u = 0.37 GeV. For the strange quarks, 
m s = 0.175 GeV, one instead has /i s = 0.46 GeV. 

Now comes the most dramatic moment, since fi u , fi s are directly connected to 
the baryon magnetic moment. Namely, similarly to the color-magnetic term, also 
interaction with usual magnetic field has the form 

5A = J2 dr k e k ^B = W 6 -^^dt, (107) 
k=1 Jo k=1 Jo 2/i fc 

where e k is the electric charge of the quark, a^ k ' is the corresponding spin operator, 
and the definition fl89|) of the constituent mass was used. 

Introducing the z-component of the magnetic moment operator 

fc = £^T' (108) 

one can write the BMM as matrix elements 

V-B = (*b\P-z\*b) (109) 



where is the eigenfunction of Q104] ), and [i k is taken at the stationary point, 
given by ( |105| ). 

For the baryon wave function we shall take here the simplest approximation, 
namely 

^ B = ^ symm (r)^ symm (a,f)iP a (color), (110) 
where ip(cr, f) is the spin- flavour part of the wave function. 
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Table 4: Magnetic moments of baryons (in nuclear magnetons) computed using 
Eqs. (|TT|) , ( |108| ) , (|i~4|) in comparison with experimental data from PDG [[24] 



Baryon 


P 


n 


A 


£- 


S° 


s+ 




H° 


Q- 


A++ 


Present work 


2.54 


-1.69 


-0.69 


-0.90 


0.80 


2.48 


-0.63 


-1.49 


-2.04 


436 


Experiment 


2.79 


-1.91 


-0.61 


-1.16 




2.46 


-0.65 


-1.25 


-2.02 


4.52 



Results of calculations of baryon magnetic moments |)5j are given in Table 4 in 
nuclear magnetons. 

Note that all our predicted values of are calculated only through the string 
tension and this is the first calculation of this sort, which tells that the string picture 
of baryons (we call it the QCD string model) is very good starting approximation 
- agreement of the calculations with experiment is on average better than of any 
other approach. 

One more example in this section will illustrate the universality of the approach, 
where the notion of constituent gluon mass enters and this quantity is explicitly 
calculated. 

The example is the glueball spectrum, calculated in this approach in [^, [7(| . 



Since gluons have no current masses, = 0, the hamiltonian for angular mo- 
mentum L = has the form (without spin and L corrections) 

P 2 

H = \-Ho + <Ta#r (111) 

Mo 

Here a a dj = \crj und , <7f und = 0.18 GeV 2 . The eigenvalues are easily obtained from 

fa(n)\ m 

M (n) = 4^ (n), /x (n) = y/oadj I — — J , (112) 

and a(L = 0, n = 0) = 2.338, a[L = 1, n = 0) = 3.36, a(L = 0, n = 1) = 4.088. 

To this Hamiltonian Hq one should add spin-dependent terms H s (totally ex- 
pressed through a s and D, DiseeW^andref s. therein) and orbital correction Hl 



derived in |58| . Referring the reader to [^9|, for details here we only give the 
Table of galueball masses in comparison with two sets of lattice calculations (for the 
same value of a/ = 0.238 GeV 2 , chosen on the lattice). 

One can notice a remarkable agreement between analytic and lattice data. Note 
also that both calculations have essentially the same input: the string tension er/, 
while spin splittings are also defined by the value of the strong coupling a s . It means 
that the Field Correlator Method is also successful in describing gluon bound system 
and the notion of valence gluon and its constituent mass (energy) is fully supported 
by independent calculations. 

Let us turn now to another gluon system, which has important theoretical rather 
than experimental meaning. This is the gluelump system calculated both on the 
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Table 5: Comparison of predicted glueball masses with lattice data (for <jf = 0.238 
GeV 2 ) 





M (GeV) 
This work 


Lattice data 


Ref. [1] 


Ref. [3] 


0++ 


1.58 


1.73±0.13 


1.74±0.05 


0++* 


2.71 


2.67±0.31 


3.14±0.10 


2++ 


2.59 


2.40±0.15 


2.47±0.08 


2 ++* 


3.73 


3.29±0.16 


3.21±0.35 


o-+ 


2.56 


2.59±0.17 


2.37±0.27 


-+* 


3.77 


3.64±0.24 




2-+ 


3.03 


3.1±0.18 


3.37±0.31 


2 -+* 


4.15 


3.89±0.23 




3++ 


3.58 


3.69±0.22 


4.3±0.34 


1— 


3.49 


3.85±0.24 




2~ 


3.71 


3.93±0.23 




3~ 


4.03 


4.13±0.29 





lattice |43| and in the framework of the present method in |44| . There is a remarkable 
agreement between the states defined in the same way (through space-like links). 

The most important feature of gluelumps is that the lowest masses with electric 
quantum numbers 1 correspond exactly to the masses observed in the exponents 
by the Pisa group |p] . And indeed calculations in Q| for the electric gluelump 



(which is described by electric correlator (E^Ek)) yield for the mass /x ~ 1 GeV in 
agreement with ||l] . Note however that magnetic gluelumps with quantum numbers 
1 H , probably could not been observed in JO]], since quantum numbers were not 



specifically projected in |41| , and therefore presumably only lowest mass could be 
observed. 

Analytic calculation of gluelumps, i.e. of field correlators with given given quan- 
tum numbers, signifies a new important period in the development of the Method 
as whole: the functions D(x) and D\(x) which have been before introduced as input 
from lattice data, will be calculated analytically within the same method, making 
finally the whole theory complete and self consistent. The first attempts in this di- 
rection have been already done in |J, where fi = 1/Tg was calculated in terms of 
string tension. 
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8 Dynamics of the temperature phase transition 
in QCD 

Calculations of the free energy of gluons and quarks, which have been done in section 
4, will serve us in this section to compute the phase transition temperature and to 
define the physical picture of the phase transition - in the first approximation of 
1/N C expansion - which will be a very realistic picture. . 

We shall consider two possible phases - the deconfined one at T > T c and the 
confined phase at T < T c . The choice of the phase from physical point of view is 
dictated by the minimal value of the free energy (or the maximal value of pressure). 
The most important ingredient in the QCD free energy (which is absent in QED) is 
the nonperturbative energy density, connected to the gluonic condensate 

= PM,p*pa\ _Hjy ^/E a E a + H a H a ). (113) 
Thus one write the free energy in the low phase (confined) as 

yF low = e - T£ ^g^e"^ - |V + 0(1/AQ (114) 

3 k 

where are massive particles (not pions) , pions are given by the third term on the 
r.h.s. of (ITU . 

For deconfined phase one should specify the condensate term ( |113j ). Since con- 
finement is electric, deconfinement would logically imply disappearence of the elec- 
tric part of the condensate (E a E a ) in ( |113| ). However the NP part of Di(x) does 
not confine and therefore in (E a E a ) = D E (0) + -Di(O) the part Df (0) could sur- 
vive. However lattice calculations of the Pisa group imply, that both D E (0) and 
Df (0) vahish near T = T c , while the magnetic part D B (0) + -Df (0) stay intact and 
contribute around one half of the total condensate, as it is for zero temperature. 

Hence one can write the free energy for the high phase as 

1 1 T 4 7T 2 7tt 2 

—F hlQh = -£- (N 2 C - 1) —N.T^rif (115) 

V 3 9 2 v c ; 45 180 i K ' 

where rif is the number of flavours, and the second and the third term refer to gluons 
and quarks computed respectively in ( |54"D and (|70|). 

Comparing (|114|) and ( |115|) one can see that Fhi g h starts higher than Fi ow at 
T = (note that e < 0) but drops faster and at some T = T c both F are equal, 
which means that the phase transition takes place and nature chooses that phase 
which ensures lower free energy (higher pressure) (this is a consequence of the 2nd 
thermodynamic law). 

This value of T r is (ITi 



U N a 8 {F) 
6 iVc 32tt 



S(^c 2 - 1) + 7 &N e n f 



1/4 

(116) 
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Note that at large N c the critical temperature is 0(N®), which is expected, but no 
other model predicts this correct dependence. 

It is instructive to compare the values of T c from ( |116| ) and from recent lattice 
calculations 
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Table 6: Values of T c for N c = 3 and different rif compared to lattice data |7I| 



SU{3) 


T c (m p = 0.77) (MeV) 


Eq.(|116|) MeV 


n f = *) 


270 


220-240 


n f = 2 


170 


150 


n f = 4 


131 


134 



*) For nf = the quenched value of the gluonic condensate was taken 3-4 times 
larger than the standard (nonquenched) value ^(F 2 ) = 0.012 GeV 4 . This in agree- 
ment with arguments of SVZ[Q] and lattice calculations of Di Giacomo et al [ |2 . 

Several comments are in order. First, values of T c are surprisingly close to 
the lattice calculations which supports our picture of phase transition. Secondly, 
magnetic condensate is as strong as for T = (and becomes stronger with grouth of 



T). This leads to the spacial correlations above T c , calculated in [2H| and measured 
on the lattice. 

Finally, one should not confuse the above physically correct picture of phase 
transition with a popular model, where all system is immerged in the huge bag 
with the standard bag pressure (which is order of magnitude smaller than gluonic 
condensate [0]). This latter picture leads to many inconsistencies and cannot be 
considered as a serious suggestion. 

Recent development of our picture demonstrating a good description of non- 
perturbative dynamics of quark gluon plasma was done for SU(2) in plj and quite 



recently for 577(3) in [72 1. Further work incorporating both electric field in Polyakov 
lines and magnetic field modifying gluon and quark free energy is now under inves- 
tigation. 
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